Estimating the scaling function of multifractal 
measures and multifractal random walks using ratios 

in 

(N 

Abstract 

h' 

^/^ ' In this paper we prove central limit theorems for bias reduced estimators of the 

(-H ■ structure function of several multifractal processes, namely mutiplicative cascades, 

multifractal random measures, multifractal random walk and multifractal fractional 
random walk as defined by Ludena (2008). Previous estimators of the structure 
functions considered in the literature were severely biased with a logarithmic rate 
of convergence, whereas the estimators considered here have a polynomial rate of 
^ \ convergence. 

\o ' 

^ ■ 1 Introduction 
O 

A random process X with stationary increments will be called multifractal if its scaling 
behaviour is characterized by a strictly concave function (, such that for a certain range 
of real numbers q 



^; E[|A(t)-X(.)n=c(g)|t-.|«''). 

If the function ( is linear, then the process is said to be monofractal, as is the case for 
instance for the fractional Brownian motion (FBM) B^, < H < 1, which is defined as a 
continuous centered Gaussian process such that Bh{0) = and 

var(5H(t) - Bnis)) = \t - s^^ . 

Then, for all q > -1, E[5H(t) - Bnis)]"] = c{q)\t - s]"^ , with c(g) = E[|5^f(l)|^] 
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Several truly multifractal processes with stationary increments have been defined. The ear- 
liest one is the multiplicative cascade introduced by Mandelbrot (1974) and rigorously stud- 
ied by Kahane and Peyriere (1976). These processes were generalized by Barral and Mandelbrot 
(2002), Muzy and Bacry (2002) and Bacry and Miizy (2003). The latter authors intro- 
duced multifractal random measures (MRM) and multifractal random walks (MRW) as 
changed time Brownian motion. Ludeiia (2008) and Abry et al. (2009) introduced multi- 
fractal (fractional) random walks which are conditionally fractional Gaussian processes. 

For all these processes, multifractality results from a distributional scaling property which 
can be written as 

{X{Xt), 0<t<T}^= {UxX{t), < t < T} , 

for < A < 1, is a positive random variable independent of the process X such 
that ^[Ul] = A^^''^ for q < g^ax a certain parameter depending on the process under 
consideration. For the models we will formally introduce in the sequel, it is defined as 

gmax = sup{g : C(g) > 1} • 

It is also important to note that the fixed time horizon T beyond which this scaling property 
need not be true is finite, except for monofractal processes such as the FBM. 

Given a multifractal process observed discretely on [0,T], it is of obvious interest to be 
able to identify the scaling function (. 

Let ti, . . . ,tN, with ti — tj„i = A = T/N be a regular partition of [0,T] (typically on a 
dyadic scale). Typically, for q < gmax, C{q) is estimated by calculating logarithms of the 
empirical structure function 

5^(X,g):=5^|AX,f , 

j=0 

where AXj = X{jA) — X{{j — 1)A). Estimators of ( can then be defined by 

, \og,{S^{X,q)) 



CN{q) : 

C^(g) := 1 + log2 



log2(A) 

SN{X,q) 



S2N{X,q) 

These estimators have been thoroughly dealt with for multiplicative cascades in Ossiander and Waymire 
(2000). The authors show that Cm(q') and Ca/(q') are consistent estimators of ({q) for q < qo, 
where qo < gmax is the largest value of q such that 

aq)-qC{q)<q + i- (i) 

For q > Cm(q') is seen to converge almost surely to a linear function of q. Moreover, 
conditional CLTs (where the limiting distribution is a mixture of normal laws) are seen to 
hold for suitably normalized versions of both estimators if 2q < q^. 
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However, the convergence rates for both these estimators are very different. The rate of 
convergence of Ca'(Q') is of order log2(A^) because of the existence of a bias term, whereas 
that of Cn{(i) is a power of which depends on (. 

In order to enlarge the domain of consistency of the estimators and obtain unconditional 
CLTS, the so-called mixed asymptotic framework has been introduced by allowing the 
number L of basic observations intervals to increase with N. In the case of multiplicative 
cascades and MRM, the processes over different intervals are independent. The observa- 
tions are X{{jL + k)A), < j < L — 1, < k < N — 1 and the estimators are now modified 
as follows 

\og2{SL,N{X,q)) 



a,^(X,g) :=! + 
a,^(X,g) := 1 + log 

with 



log2(A) 

SL,NiX,q) 
SL,NiX,q) 



L-l N-1 



L+k 
j=0 k=0 

The mixed asymptotic framework for multiplicative cascades has been recently developed 
in (Bacry et al. (2010)). The authors show that if L = [2"^^], where [x] stands for greatest 
integer m < x with x > 0, then Cn,l{X, q) is consistent for q < q^ where q^ is the largest 
value of q such that 

C{q)-qC{q)<q + x + i. (2) 

Note that as x tends to infinity, might be greater than gmax, so we will only consider 
values of x such that < gmax- 

However, once again, there exists a bias term bn '■= W,[Mf]/n, which entails slow convergence 
of the estimator. In analogy to the non mixed asymptotic framework it is reasonable to 
consider ratio based estimators such as (N^L{X,q) in order to improve convergence rates. 
It turns out, as follows quite straightforwardly from the results of Bacry ct al. (2010), that 
CN,L{X,q) — )■ C{q)y a.s. for a dyadic partition, but the authors failed to prove a CLT 
(for 2q < q^). Almost sure convergence for dyadic partitions, or in probability for general 
partitions, of ^]\f,L{X, q) has also been recently considered by Duvernet (2009) for x ^ and 
X a Brownian MRW or a MRM. However the author does not prove CLTs nor establish 
convergence rates in either case. 

The main goal of this paper is to obtain CLTS for the estimator C, in the mixed asymptotic 
setting, for multiplicative cascades, multifractal random measures (MRM) and multifractal 
random walks (MRW) for H > 1/2. Our main results in all these cases state unconditional 
CLTS with polynomial rates of convergence, contrary to ({q) which can only achieve log- 
arithmic rates of convergence, and to the case L = 1 where only conditional CLTs can be 
obtained. 
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We will consider multplicative cascades in Section 2, MRM in Section 3, and MRW in 
Section 4. Section 5 contains the technical parts of the proofs. To the best of our knowledge 
our results are the first to deal with the MRW in the case H > 1/2. 

Properties of log-Laplace transforms 

We conclude this introduction by gathering certain convexity properties of Laplace trans- 
forms that we will need. In all the models we consider, the function C, can be expressed 
C{(1) = (1 ~ "^{q) for ^ function ip which is the log-Laplace transform of some random 
variable Y, i.e. 

^(g) = logE[e^^], 

with 1^(0) = = 1. The function tp is thus strictly convex and we can express go and 
gmax in terms of ip- 

go = max{g : q^'{q) - ip{q) < 1} , g^ax = max{g : ip{q) < g - 1} . 

The convexity of ip and = implies that g^ax > 1 if and only if ip'il) < 1, and 
■^'(O'max) > 1- This in turn implies that 1 < go < gmax- Indeed, if g > gmax, then ip'{q) > 1; 
the strict convexity of ip and = implies that ^^(g) < ip'{q){q — 1), hence 

thus g > go. Also, < 1 implies that go > 1, since ip'^l) — ip{l) = ^/''(l) < 1. 

Let now q^^ be defined as the largest g such that qip'{q) —ip{q) < 1 + X- Since ip is convex, 
the function q'ip\q) — ip{q) is increasing, thus q^ > go. 

For g > 1, we will also be interested in the positive and increasing function p i— )■ ip{pq) — 
pip(q). If pq < q^, then, by convexity, 

< ip{pq) — pip{q) = vPipq) — V^(g) — {p — < (j> — ^)pq'^'{pq) — {p — ^)i^{pq) 

= {P- - ^{pq)} < (|> - + x) • (3) 

For p = 2, if 2g < g^, this yields 

0<ij{2q)-2i;{q)<l + x- (4) 

2 Multiplicative cascades 

In this section we give a precise formulation of consistency results for C{q), whenever 
q < qxi and a CLT whenever 2g < g^, in the case of multiplicative cascades. The 
results are a straightforward application of previous results of Bacry et al. (2010) and 
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Ossiander and Waymire (2000). However, they provide the basic framework for deahng 
with both MRM and MRW so will be dealt with in some detail. 

Before we state the main results we shall introduce the basic framework for mixed asymp- 
totics following Bacry et al. (2010). For any given n-uplet r and i < n set r\i = (ri, . . . , r^) 
and if s is an i— uplet and v an n — -i— uplet set r = s * w to be the resulting n— uplet 
obtained by concatenation. 

For each j E and fixed T, set := [j'T, (j + 1)T]. Over each I^^^ we will construct an 
independent multiplicative cascade as defined in Mandelbrot (1974). For this consider a 
collection {Wr'^\r G {0, l}"r2 > 1, j G Z} of independent random variables with common 
law W such that E[W] = 1 and E[W^log2 W] < 1. For each n > 1 and j E Z, consider the 
random measure defined by 

n 

A(f)(/)=T2-" rw^' 

rG{0,l}"n/ j=l 

for any / a Borel subset of over and each r = (ri,...,r„) E {0,1}"" is iden- 
tified to the real number *"i2"~^'. It can be seen (see Kahanc and Pcyricrc (1976), 
Ossiander and Waymire (2000) for details on the construction and main results) that there 
exists a random measure such that 

P(A(/') ^ A^;^,) as n ^ oo) = 1 , 

where =^ stands for vague convergence. The limiting measure verifies E[A^^([0, T])] = T. 
By construction A^'' are independent random measures, defined over the disjoint intervals 
/(^■). Set Aoo:=E,gzAS. 

Set = a{Wr^'\r E {0, 1}", j E Z} and let A^f'^ := [(j + k2-^)T, (j + {k + 1)2-")T], 
A; = 0, . . . , 2" — 1 be the k—th. diadic interval at level n, of the interval Then, 

n 

Aoo(Al^l)=2-%,,,„ U^i^H, (5) 

where for each n, Zj^^^n, < k < 2", j E Z, are i.i.d. random variables with the same 
distribution as Aoo([0, 1]) and independent of and r„(A;) is the dyadic representation of 
k, i.e. k = X]r=i ^ ^ 2"- Moreover, Zj^2k,n+i and Zj^2k+i,n+i are independent 
of Zj^k',n for k' 7^ k. The identity (5) straightforwardly yields the scaling property: 

E[AL(Ai-';i)]=2-"'^^^^E[A^^([0,l])]. 

with 

aq) = q-\og2{E[W^]) . 

It is shown in Kahane and Peyriere (1976) that for q > 1, the condition ({q) > 1 implies 
E[A^^([0,T])]<oo. 
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Example 2.1. Consider the log-normal cascade, where log2 W = fi + aZ and Z ~ N(0, 1). 
The condition K[W] = 1 implies that fi = — cr^/2. Then it is easily obtained that 



^(g) = g(g - l)aV2 , g^ax = (2/a2) V 1 , qo = V2/a, = ,/{2/^W+)d ■ 
Denote 



j=0 k=0 

and 



C(g) 



, log2 G^L.n g . ( SL,n{q) \ 

n \SL,n+l{q)J 



Consistency 



For each n > 1, let ^j,k,n, 0<j<L — 1,0<A;<2" — l}bea collection of i.i.d. random 
variables, independent of Define 

L-l 2"-l n 
j=0 A:=0 i=l 

In Bacry et al. (2010) the following general result is shown to hold. 

Proposition 2.1. Assume that q < qx there exists e > such that E[^^+''] < oo. If ^ 
is non negative, then 

L-'2-"2<('')(^„,g-E[5„,,])^0, a. s. 

Note that by construction E[5^„^g] = L2"2~'^''*^'^)E[,^], so that the above result yields the 
convergence L~^2~^2'^^^''>Sn,q — )■ E[(^], a.s. under the stated conditions. As a consequence, 
by the definition of Si^niq), Proposition 2.1 yields 

^-i2-2<('')^z.,„(g) ^ E[A^([0, 1])] (7) 

a.s. for q < q^- Then, clearly, 

n 

and this implies that ({q) — )■ ({q) a.s., so that 

C(g) ^ C(g) a.s. 

On the other hand, if g > q^, then Bacry et al. (2010) show ({q) — )■ C{qx)q, which is a 
linear function of q. In this case ({q) is also not consistent as the normalized structure 
function tends to zero (Ossiander and Waymire (2000)). 
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Central limit theorem 

Based on Proposition 2.1, it is also possible to obtain a CLT for C{q). We remark that 
in the mixed asymptotic framework the limiting variance is deterministic. The proof of 
the CLT follows from a series of corollaries of the following general result for the mixed 
framework which is a direct generalization of Proposition 4.1 in Ossiander and Waymire 
(2000) and Proposition 2.1. We first state some general notation. Let {^,^^,^,^0 < j < 
L — 1,0 < k < 2"~^, 7T, > 0} be as above and let Sn,q be as in (6). Define now 

j=0 k=0 i=l 
Rn,q ■— ^n,q/'^n,q ■ 

The following proposition is seen to hold true as a direct generalization of Proposition 4.1 
in Ossiander and Waymire (2000), whenever 2q < q^. 

Proposition 2.2. //E[e,-fc,n] = 0, = and if 

sup sup E [\ij,k,n?^^^^^] < , 
n j,k 

for some 6 > 0, then for 2q < q^ under the assumption M([0, 1]) > 0, a.s. 

lim E [e*"^-' I J-J = e""'^'/' 

n— 5>oo 

and Rn^q converges weakly to the centered Gaussian law with variance a^. 

The proof follows exactly as that of Proposition 4.1 in Ossiander and Waymire (2000), 
using Proposition 2.1. The latter also yields that L~^2~^2'^^^'^'^'>Vn^q converges to 1 a.s. We 
now have 

Proposition 2.3. If 2q < q^, then 

^-i/22-n/22nC{2,)/2 _ 2^^'')-' S L,n+l{q)] N{0, V {q)) , 

with 

V{q) = var [Zl - 2^^'?)-^-'? {ZlWf + ZlW^]) 

and Zl, Z2 are i.i.d. with the same distribution as Aoo([0,l]) and independent ofWi, W2 
which are i.i.d. with the same distribution as W and Zq = {ZiWi + Z2W2)/2 has the same 
distribution as Aoo([0, 1]). 



7 



Proof. The proof follows from Proposition 2.2, by noting that SL,niq) — '2'^^''^ ^SL,n+iiq) 
can be expressed as 



j=0 k=0 i=l 

with 

c — yi _ nC(q)-l-q / yl i 79 1^/9 \ 

'ij,k,n — ^j,k,n ^ |_^i,2fc,n+l*^''^-rn(fc)*0 ^ j,2k,n+l^'' rn{k)*l j 

since r„(A;) * = r„+i(2A;) and r„(fc) * 1 = r„+i(2A; + 1). Indeed, the r.v.s C,j^k,n, j ^ ^5 
< A; < 2" are i.i.d. (for each fixed n) and it clearly holds that IE[0,fc,n] = 0, IE[^|,fc,n] = ^i^l) 
and E[|^j^fc^„p+^] < 00, whenever 2q < q^ax for small enough 5 > 0. □ 

Thus we obtain. 

Theorem 2.4. Assume 2q < q^. Then 

2n(i+x+2^(,)-^(2.))/2|^~(^) - aq)} iV(0, r(g)/E[AL([0, 1])]). 

Proof. By Proposition 2.1 and (7), S^n+M'^^^^"'^"^ / -> 1 a.s. so 

C(g) - Cl?) = iog2 ^...wi o 73 = - iog2 1 - 



^i,4g)-2^('?)-i5i,„+i(g) 



+ {l + op(l)} . 



The proof is concluded by applying Proposition 2.3 and noting that 2"*^^+^"*"^*^'^'' '^(2''))/2 ^ 
L2"2~"'^('^)/L^/^2"/^2~"^(^^)/2. □ 



3 Multifractal random measure 

Once again we are interested in the mixed asymptotic framework defined by the param- 
eter X- The main ideas dealt with in this section are very similar in spirit to those in 
Duvernet (2009), we include the proofs for completeness' sake, since they are very similar 
to those which will be developed to study Multifractal Random Walks. 

We recall the main definition and properties of Multifractal Random Measures, hereafter 
MRM, foUwing Bacry and Muzy (2003). 

Start by defining for / > 0, wi{u) = P{Ai{u)) and set 

Mil) = lim [ e'"'^") du , 
1^0 J J 
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where / is any Borel set in M. Here P is a set process over = {{s,t),t > 0} such that 
P{A UB) = P{A) + P{B) and P{A) and P{B) are independent if A n = 0, and 

E[e'?^(^)] = e^^'')'^^^) (8) 

with fi{A) = J^t~^ ds dt and 

Ai{u) = {{s,t),u- {t/2 AT/2) < s <u + {t/2AT/2),t> 1} . 
It will be useful to note that 

fi{Ai{t)) = T + \og{T/l) . 

The function ip is the log-Laplace transform, assumed to exist for q < q*, for some q* > 1, 
of the infinitely divisible random measure P. It is convex and satisfy ip{0) = = 0. By 
the Levy Khinchine representation Theorem, it can be expressed as 

^2 roc 

^(9) = — + mq+ / {e'^'' - 1 - xl{|^|<i}}z/(da;) , 
^ ^-00 

where u is the Levy measure of P and satisfies 




(x^ A l)z/(dx) < 00 . 



The assumption that ipiq) is finite for q < q* entails the following condition. For all q < q*, 

/oo 
e'^^z/(dx) < 00 . 



Aiiu) 




Figure 1: The set Ai{u) 
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By Theorem 4 in Bacry and Muzy (2003), there exists a certain infinitely divisible r.v. Q\ 
which is independent of M[0,T], such that E[e'?^^] = A^^^*) and for X,l e (0, 1), 



{wxi{\u) , < n < T} = {wi{u) + Qx, < u < T} . 



(9) 



The latter is known as the scaling property. This implies that 



M([0,AT]) = Xe^'M[0,T] 



(10) 



for A G [0, 1], so that 



E[M'?([0,AT])] = A'^('?)m(g) 



(11) 



with ({q) = q — ip{q) and 



m{q) = E[M''[0,T]] . 



It can be seen (Bacry and Muzy (2003)) that C(g) > 1 ^ E[M'J{[0,T])] < oo. As previ- 
ously, set gmax to be the greatest value of q such that C(g) > 1 and for x ^ 0, define q^ as 
the greatest value of q for which qip\q) < ip{q) + \ + x- Assume moreover that x is such 
that qy^ < gmax and recall that by convexity, if g < g^^, then 



Example 3.1. Consider the Poisson cascade introduced by Barral and Mandelbrot (2002). 
Let be a Poisson point process with intensity measure /i on (— oo, oo) x (0, oo]. Let Pj, 
z G Z denote the points of A^ and let {IV, PVi} be a collection of i.i.d. positive random 
variables such that W\W] = 1. Define the random measure P by 



for all relatively compact Borel sets A G (—00,00) x (0, 00]. Then (8) holds with ip{q) = 
ElW^] - 1. 

g„,ax = max{g : E\W^\ < g} , q^ = max{g : gE[lV'^(log(lV) - 1)] < 1 + x} • 

Example 3.2. The random measure P can be a Gaussian random measure. Then P(A) ~ 
N(— (T^/i(A)/2, cr^/i(/l)) and '^{q) = cr^q{q — l)/2 so that we get the same values of gmax, 
go and g^ as for the multiplicative cascade of the previous section. Note that in this case, 
var(P(y4)) = ^"(0)/i(A) is finite if and only if n{A) < 00. 

Example 3.3. Let a G (0, 1) and P be a totally skewed to the left a-stable random measure, 
i.e. -{/^(g) = 0""(g — g"). Then g^ax > 1 if and only if (t"(1 — a) < 1 and then g^ax = 00 
and for x ^ 0, g^^ = cr~^((l + x)/(l ~ a))^^". It is noteworthy that contrary to the 
previous case, we have here that E[|P(74)|] = 00 for all A such that fi{A) > 0, though 
E[\P{A)\P] = Cp^^aP^i{Af/'^ i{p<a and ^i{A) < 00. 



^(g)-2V;(g/2)<l + x. 



P(A) = 5^1og(H^,)l{r,eA} 
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Example 3.4. Let a G (1, 2) and P be a totally skewed to the left ct-stable random measure, 
i.e. ipil) = — q). Then g^ax > 1 if and only if (T"(a — 1) < 1 and then g^ax < oo. 

For X > 0, = a~\{l + - 1))'/". 

With this notation, define as in the previous section L = [2"^], A^"']^ = [{j + A;2~")T, (j + 
(A; + 1)2-'^)T] and 

j=0 k=0 



CM{q) = 1 + logs 



SL,n+l{M,q) 



Consistency 

Fro convenience, denote r(g) = 1 + C{q)- 
Proposition 3.1. For q < q^, 

L-^2^-('^)SLAM,q)^m{q) , a.s. 

Plugging this into the definition of Ca/(q') yields the consistency of (Aiiq)- 
Corollary 3.2. For q < q^, 

Chdiq) C{q) , a.s. 

Central Limit Theorem 

We next give a central limit Theorem for Cm(q') in the mixed asymptotic framework. Define 
the centered random variables 

D,,,,.,, := M'^(Agi) - 2-('^)(M'^(A(i?,„^,) + MnA(^^,,„^,)) (12) 

and Dj^n,q = /_^k=o ^j,k,n,q- By construction, the variables Dj^k,n,q are centered, and we 
can also bound their covariances. By stationarity and 2-dependence with respect to j, we 
can consider the case j = 0. The following bounds for the covariances is proved at the end 
of section 5.1. 

Lemma 3.3. // 2g < gmax? then for k = 1, . . . , 2"-~^, 

2"«2.)E[Z}o^o,n,,/^o,.,n,,] = 0(A;-{^(2^)"2'^(^)+i>) . (13) 
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We will start by proving a CLT for 



Since the random variables Djn^q, < j < L — 1 are 2 dependent, it suffices to show that 
for some p > 1, 

lim = . 14 

Although the variables Dj^k^n,q are conditionally i.i.d., unlike multiplicative cascades, 
they are not conditionally centered. Hence it is not possible to repeat the proof of the CLT 
for multiplicative cascades and it follows that it is necessary to require that x > 0- 

We ffist need an expansion of E[Dq^^^^]. Set 



dq = E 



(M'^([0,T] -2^(«){M''([0,T/2]) + M'?([r/2,T])})' 



and 4,, = 2<(2g)E[Doo„_gDo,„ J. By the scaling property, E[Dl^^^^^] = 2-<(^'^^dq and 
dk^q does not depend on n. Then 

2"-l 

E[Dl^J = 2-"^(2,)^^ ^ 2 ■ 2-"^(2«) 5^ (1 - k2~")dk,q . 

k=l 

Since ip{2q) — 2ip{q) > 0, Lemma 3.3 implies that the series Yl is convergent, so the 
Cesaro mean above has a finite limit. Thus, there exists a constant Qg such that 

lim 2-(2'')E[Do%,J = 0, . 

Next we prove (14) for p = 2, i.e. we compute the fourth moment of DQ^n,q, which exists if 
Aq < q^. By computations similar to those that yield Lemma 3.3, we can prove that 

The above discussion leads to the following result. 
Proposition 3.4. //4g < g^, then 

Dn,q = — , >d iV(0, 1) , 



or equivalently, 

L~y^2^ri2q)/2jj^^^^^ iV(0, 6,) . 
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We can now prove the asymptotic normahty of tm{(i) and Cm(9)- Denote 

_ L-^2"-(^){^z.,n(M,g) -2-('?)^z.,n+i(M,g)} 

By Proposition 3.1, Rn = o(l), almost surely, so a second order Taylor expansion yields 

' gL,n(M,g) \ 

2ril)SL,n+l{M,q)J 

, A , L^^2-('^){g^,.(M,g) - 2-('?)g^,„+i(M,g)} ^ 



CM(g) - C(g) = logs 



L-i2"-(^){^i,„(M,g)-2-('')^z.,n+i(M,g)} , ^ ^ 



L-i2-W^i,„(M,g) 

Note now that ^^^^(M, g) - 2^('')S'l,„+i(M, g) = Ylj=i Yl'j=o^ Dj,k,n,q, so Proposition 3.4 
yields the next result. 

Theorem 3.5. //4g < then 



m(g) 



4 Multifractal random walk 



Throughout this section, the MRM M and the process {wi{u)} will be as defined in the 
previous section. A multifractal random walk (MRW) is the process X obtained as the 
limit as / — )■ of the integral /q e"''^"-* dBniu) where Bh is a standard fractional Brownian 
motion independent of M; see Abry ct al. (2009); Bacry et al. (2001); Bacry and Muzy 
(2003); Ludena (2008). Recall that Bh is a continuous Gaussian, centered process with 
Bh{0) = and 

vaT{BH{t) - Bnis)) = \t - si'"" , 

for all t, s G [0, 1]. For H = 1/2, B1/2 is the standard Brownian motion and will be simply 
denoted by B. This means that X is the conditionally (with respect to M) Gaussian 
process whose covariance function is defined in (15) or (16) below according to whether 
the Hurst parameter of the fBm H = l/2orif> 1/2. Except for the case H = 1/2, which 
is ordinary Brownian motion, it is worthwhile to remark that this conditionally Gaussian 
process X is not the time changed process BH{M[0,t]). 

Throughout this section -^m will stand for conditional convergence in distribution given 
M and E^^ and varjv/ stand for the conditional expectation and variance given M. We 
consider the following two cases. 
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Case H = 1/2 Bacry et al. (2001); Bacry and Muzy (2003). The MRW X is defined 
as the centered, conditionally Gaussian process with conditional covariance 

r(s,t)=lim / e'^'^^Mu = M(s At). (15) 

The scaling function is Ci/2(q') = C('?/2), since by (10) and (11), for A G (0, 1), 
{X(At), < t < T} '= \^/^e^^'^{X{t), < t < T} , 

nnm = mMum = c.^M'/^it)] = c,m(g/2)t«^/2) 

where = E[|N(0, l)]-?] and m(g) = E[M^([0, 1])]. 

Case if > 1/2 Abry et al. (2009); Ludcna (2008); Muzy and Bacry (2002). The 
MRW X is defined as the centered, conditionally Gaussian process with conditional 
covariance 

/•* r e-'We-'(-) ^ f M(dn)M(dtO 

where Ch = H{2H — 1). This process is well defined whenever H — ip{2)/2 > 1/2, cf. 
Ludeha (2008). Convexity of ip yields ip{2) > 0. The scaling function (h is defined by 

Cniq) = qH - ip{q) , 

since by (16) and (10) we have 

{X{Xt), < t < T} '= X^e^^{X{t), < t < T} , 



with 



T pT ^ 9/2" 

\u-v\^"-^M{du)M{dv] 

Jo 



(17) 



Since we are considering the mixed asymptotic framework, we assume we have a collection 
of MRM M^^\ j = 0, . . . , L — 1, which are independent, defined over consecutive intervals 
of length T. For j = 0, . . . , L — 1 and k = 0, . . . , 2""-*^, define AXj^k,n = ^{j+(k+i)2-")T — 
X{j^k2-~^)T- As above, we will investigate the asymptotic properties of fx{q) defined by 



Tx{q) = log2 

where now 



SL,n{X,q) 
Sl ,n+l 



SL,n{X,q)=}^}_^\AXj,k,r. 
j=0 fc=0 
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It will appear that Tx{q) is an estimator of Tniq) defined for i7 > 1/2 by 

rniq) = Cniq) - 1 • 
Thus we define an estimator Cx{q) of the scaling function (niq) by 

Cx(g) = 1 + fxiq) = 1 + log, f /"'"^f;.'\ 

Denote T„(X, q) = 2^^ ('^'^ Sl, n+i{X, q) - Sl^X, q). Then 

Tn{X,q) 

SL,n{X,q) 



fx{q) -Tniq) = log 1 



We will prove that Tn{X,q)/SL,n{X,q) — 0, a.s. so that a Taylor expansion is valid and 
yields 

fx{q)-rn{q) = -^^^{l + o{l)). 

JL,n{X,q) 

In order to study the ratio above, we will first prove that if if = 1/2, then 

L-'T^^'^/^^SL,n{X,q)^Cgm{q/2) 

and if ii > 1 /2 then 

L-'r-«('')SL^^{X,q)^CgmH{q) 

with mniq) as in (17) and Cq = E[|N(0, 1)|''] in both cases. To study Tn{X, q), we write 

T„(X, q) = Tn{X, q) - Em[T„(X, q)] + Em[T„(X, q)] . 

We will prove that in both cases, T„(X, g) — ¥.M\Tn{X,q)] and ¥,M[rn{,X,q)] converge 
jointly to independent centered Gaussian distributions with the same normalization. This 
will yield the asymptotic normality of Cx{.q) — Cniq)- 

Because of the different nature of the conditional dependence structure, which yields dif- 
ferent scaling functions we will consider the cases H = 1/2 and H > 1/2 separately. 

4.1 The case // = 1/2 

In this case, it holds that 

EM[SL,n{X,q)] = CqSLAM,q/2) , 
vaiM{SL,n{X,q)) = alSL,n{M,q) . 
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By Proposition (3.1), if g < q^, we get 

L-^2''-(''mEM[SLAX,q)] -> c,m(g/2) , a. s. 
L-i2"^(^)varM(^L,n(X,g)) ^ ajm(g) , a.s. 

This implies that L~^2''^^i/^^ SL,n{X, q) converges in probabihty to Cqm{q/2). Since Si^niX, q) 
is the sum of L2^ conditionally independent terms, by Borel-Cantelli arguments similar to 
those used previously, almost sure convergence also holds, i.e. 

L-^Triil^'>SL,n{X, q) ^ c,m(g/2) , a.s. (18) 

Using the notation (12) of the previous section, we have 



EM[T„(X,g)] = Cq2<''/^^SL^n+i{M,q/2) - c,5z.,„(M, g/2) = -c, J] ^ 

j=0 fc=0 



L-1 2"-l 

n. . 

n,q ■ 



Thus, by Proposition 3.4, if g < then L-i/22"^(9)/2Ea^[T„(X, g)] converges to a centered 
Gaussian random variable with variance S(l/2,g) = 0^0^/2- By the conditional inde- 
pendence of B and M, T„(X, g) — Em[T„(X, g)] is a sum of centered and conditionally 
independent random variables with conditional variance 

varM(T„(X, g)) = aj5^,„(M, g) + ^^^(a^-C'?/^) _ 2-(''/2)+i)5z.,„+i(M, g) . 

where = var(|N(0, 1)|^). By Proposition 3.1, L~^2'^'^^'^\aj:M{Tn{X^q)) converges almost 
surely to the positive constant r(l/2,g) defined by 

r(l/2,g) = ajm(g) {l + {2''<'^/^^ - 2-to/2)+i)2--(^)} . 

Thus, L-^l'^2''<'i^/^{Tn{X,q) - Em[T„(X, g)]} converges weakly conditionally on M to a 
Gaussian random variable with variance r(l/2,g), independent of M. Since the variance 
is deterministic, this assures non conditional convergence to the stated Gaussian r.v. More- 
over, the conditional independence of B and M also implies that the sequence of random 
vectors 

L-y^2^-i^)/^{Tr,{X, q) - Em[T„(X, g)], Em[T„(X, g)]) 

converges weakly to (Zi, Z2) where Zi and Z2 are independent centered Gaussian random 
variables with respective variances r(l/2,g) and S(l/2,g). 

The previous considerations yield 
Theorem 4.1. If q < q^, then 

r(l/2,g) + S(l/2,g) 



^l/22n(^(,/2)-^(,)/2+l/2)^^^(^) _ iV 



c>2(g/2) 
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4.2 Case F> 1/2 

We must first obtain an equivalent of (18). Tliis is done in the following Lemma whose 
proof is postponed to Section 5.2. 

Lemma 4.2. 

• Ifq < qx ihen L-i2"^«('')EM[^L,n(^,g)] Cqmniq) a.s. 

• If2q< and H < 3/4, there exists a positive constant T{q) such that 

i^-'22"-«(«)varM(5i,n(X,g)) ^ r(g) , a. s. 

Moreover 

2,-i2-H(.)5^^^(X,g) ^ Cgmniq) , a. s. (19) 

Consider now the conditional variances 

aj,k,n,H = E]^^[{AX^,k,n)% for H > 1/2 
and the conditionally standard Gaussian random variables 

Let Gg{x) = \x\'^ — Cg. With this notation 
T„(X,g)-EM[T„(X,g)] 

L-l 2"-l 
j=0 k=0 

Once again, arguing as in the proof of Lemma 4.2, one has the convergence 

varM(2"(2'^('')-'^(2'^)+^+^)/2{T„(X,g) -EM[T.(X,g)]}) ^ T{H,q) , a. s. 

with r{H,q) a certain positive constant. Following the proof of Theorem 3.1 in Ludcha 
(2008), the latter together with the bounds for the joint covariance structure of the con- 
ditional gaussian random variables Yj^^^n establish the following result. Ifl/2<iJ<3/4 
and 2q < q^, 

^-i/22-n.«(2.)/2|^^(^^ g) - EmIUX, q)]} N{0, T{H, q)) (20) 

As for the case H = 1/2, the fact that T{H, q) is deterministic establishes non conditional 
convergence in distribution. 
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Let us now study Ea/[T„(X, g)]. We have 

Em[I;(X, q)] =C,J2Y1 {'^^"^'^{<2k,n+l,H + 42fc+l,n+l,//} " ^Uu^h) ■ 

Denote Un,k = 2^«(''Hao,2fc,n+i,H + «o,2fc+i,n+i,//} " «o,fc,n,H- Then ?7„,fc is centered and 
'd{q) = 2"''^*^^'^^var([/„ o) does not depend on n. Similarly to Lemma 3.3, the following 
covariance bound holds. 

Lemma 4.3. 

2-Ch(29) |^[^^^^^ jj^^^^ I < ^^-{V-(2g)-2V.(g)+l} _ ^21) 

We can now compute var(^^^g"'^ Un,k)- By stationarity, 

(2"-l \ 2"-l 

k=0 / k=l 

The series cov{Un,o, Un,k) is convergent, thus the Cesaro mean above converges to its sum, 
and we obtain that there exists a constant q) such that 

Iv"'2""^«(2'?)var(EM[T„(X,g)]) ^ , 

This and (20) yield the asymptotic normality of the estimator. 
Theorem 4.4. If Aq < q^ and < 3/4 then 

2n(i+x-^(2,)+2,(,))/2^^^^^) - Cniq)} N fo, nH,q) + m,q) 

5 Proofs 

In all the proofs, without loss of generality, we set T = L We preface the proof by stating 
some results for infinitely divisible random measures. The infinitely divisible measure P 
introduced in Section 3 can be decomposed as P = Pq+Pi where Pq and Pi are independent 
and 



with 

^/'o(g) 



^2 roo 

—q^ + mq + j {e'^'^ - I - qxl{\^\<i}}u{dx) 

I {e''" - l}z/(dx) . 
J —00 
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Then, for A such that fi{A) < 1 and g > 1, it holds that 

E[|e^i(^) = 0(/i(A)) , (22) 

E[\e''°^^^-l-Po{A)\'^] = 0{fi{A)). (23) 

Further, write 

e^(^) _ 1 - Po(A) = {e^i(^) - l}e^o(^) + e^^^^) - 1 - PoiA) . (24) 
This decomposition, (22), (23) and the independence of Pq and Pi yield 

E[|e^(^)-l-Po(A)n = 0(/i(^)). (25) 

Since P, Pq and Pi are independently scattered, these inequalities yield martingale maximal 
inequalities. 

Lemma 5.1. For A such that fi{A) < 1, and for Cu an increasing sequence of measurable 



subsets of A, it holds that 

E[supjPo(a)|] = 0(//2(A)), (26) 

E [sup, le^^^'^) - 11"] = 0(/i(A)) , q>2, (27) 

E [sup, |e^(^") - 1 - Po(a)r] = 0(/i(A)) ,g > 1 . (28) 



5.1 Multifractal random measure 



Proof of Proposition 3.1. For e > such that (1 + e)q < q^ set e' = min(2£:, 1/2, x), 
no = [1/e'], a = 1/no and finally In = 2-(^-")". Let 



L-l 2"-l 



We will prove that there exist some constants C, 77 > 0, such that we have 



E 



T 

-J- 7 



1 



< C2 



-nr) 



E 



T 



^L,n(M,g) 



E[r„J E[5i,n(M,g)] 



< C2' 



-nr) 



(29) 
(30) 



The above inequalities and an application of the Borel-Cantelli lemma yield that 



T 



n,q 



and 



E[T„,,] 
SlAM, q) 



— )■ 1, a.s. 



T 



n,q 



— !■ 0, a.s. 



E[SL,n{M,q)] E[T„,,] 
Since E[^L,„(M,g) = L2-"^(«)m(g), Proposition 3.1 follows. 



□ 
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Proof of (29). Set < x' ^ X ^ fo^' as defined above and let 5 > be such that 
(1 + 5)q < q^i. Now set e = inf (5, 1/4, x/2). Since x' < convexity of the function ip 
assures e < 5 < e and hence by construction we have that 2e < a. 

We have, for all j, k, n, 

E[M5(Ag;,)] = 2-"«^)m(g) , 

so that E[SL,niM,q)] = L2-"^(«)m(g). Next, by (8), 

E[T„,q] = L2"E[e'''"'»(^'+'^2"")] = L2''e^^''U-^^'''^ . 

The variables e^"''"*^-'"'"^ — E[e '•■'"'"^ "'^•'j are 2— dependent (in j) and centered, so there 
exists a constant C > such that 



E 



T 



n,q 



< — E 



2" — 1 



On 



fc=0 



-'/'(g) 



- 1 



Applying Lemma 5.2,we have that for any e' < e, that there exists a constant C such that 



E 



T 



1 



< ^2"[-^X+(l-"){'/'((l+e)9)-(l+e)^(9)-f'}] 



By (3), for all e > such that g(l + e) < g^, we have 

- ex' + (1 - a)mi + e)g) - (1 + e)^(g) - e') 

< e[-x' + (1 - a)(l + x')] - (1 - = -"^(1 + x') + (1 - «)(e - e') . 

This can be made negative since since e' can be chosen arbitrarily close to e. This 
proves (29). □ 

Proof of (30). We start by using again the argument of 2-dependence in j, to obtain, for 
e and for some constant C, 



E 



1 ^pV- 

Z2" ^ ^ 2^ 



j=0 k-- 



< -E 



2"-l 



on 



fc=0 



For clarity, we now omit the superscript (0) in A[,°^. Let M„ denote the random measure 
with density e"""'" with respect to M . By construction, the measure M„ is independent of 
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the process wi^ and 



{e 



l}M„(diz) . 



This yields 



+ 



2""C('?)m(g) 



2-<(9)m(g) 



By (39) and (40) in Lemma 5.3 below, we obtain 

E[M„^(Afc,„)] = e-'^(^)2-'^«(«)/f '?)m„(g) 
with lim„_^oo ''^n(Q') = rn{q) and 

|m„(g)-m(g)| <C2-"V 



Thus we can express Vn,k as 



2-<('?)m„(g)e-'/'('')/f^) 



mn{q) 



"^n(g) 

m{q) 



_ ^ I "^n(g) ^ y»n(g) _ ^ 

m(g) m(g) 

,gu,,„(fc2-") 



,qt«,„(fc2-") 



-V'(g) 



- 1 



Denote p = 1 + e. By Minkowsky's inequality, we get 



2"-l 



-E 



A/5 (A 



(0)^ 
k,n) 



k=0 



l+e' 



< 



+ 



W2n(g) 

m{q) 



2"-l 



on 



fc=0 



"^n(g) _ ^ 
m{q) 



^ 2"-l 



ggt«i„{fc2-") 



-E 

On / ^ 



fc=0 
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The random variables M„(Afc,n) are 2"7„ dependent and e""'" is independent of M„, thus 
we have 



E 



2"-l 



4 E 



fc=0 



< -T, — 



E 



qPQWi^ (0) 



- 2(1-Q)<:" 

The bound (3) yields ip{pq) — vi'i.l) ^ ^ (1 + x'); thus 

ilj{'pq) - pi){q) - (1 - a)e - ex < e(x' - X + ") < 

because of the choice of x' < X ~ We now bound the term in (33) by applying the 
bounds (31) and (38). We have 

2-1 (fc2-") 



rriniq) 



m(q) 



-E 

on / ^ 



k=0 



For the sum in (34), we use Jensen's inequality and the bound (41) of Lemma 5.3 to obtain 



E 



2"-l 



fc=0 



Pi 



< CE 



2-nC(g) 



^ 2-na/22n{ip{pq)-pip{q)} 



(36) 

In order to prove that the bounds (35) and (36) are good enough, we must now check that 
i^ipq) ~ VPiq) — as — ex < for small enough e and p = 1 + e. Indeed we have proved 
earlier that il^ipq) — inp{q) < e(l + x); so 

'^{pq) ~ P'^iq) — — ex < e — a/2 < , 

since e < a/2. □ 

Lemma 5.2. Let a = I/uq for some arbitrary integer uq > 2. For all p > 1 such that 
]g[gP9«'i(o)j ^ g^jiy g (0,p — 1), there exists a constant C such that 

"1 e'^"''"*^") 



E 



■ du 



E 



2-5: 



2 ^qwijk2-") 



< (7/~{'^(P9)-P'/'('?)-e'} 

< (7/~{'^(P'?)-pi/'(g)-<;'} 



(37) 
(38) 



Proof. The choice of a implies that (1 — a)nQ = rio — 1 is an integer. Denote gn{u) 
eij'^i„(«)/]E[gij'"i„(o)]_ pix some integer ko, and define rii = koUQ. If rii < n, then 



gn{u)du= / gn,{u)du+ / {^(^(n) - 5(„j(u)}du 





2(1-")"! —1 



/ gn,{u)du+ / {gn 



U) - gni{u)}du 
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We bound the first integral by applying Jensen's inequality: 



E 



(u) du 

Since wi„_^ is independent of wi„ — wi^_^ , we can write 

Thus we see that the integrals /. {fi'n(w) —gm (w) } du are centered and 2-dependent condi- 
tionally on J-'^ the sigma- field generated by {wi^^ {u),u G [0, 1]}. Thus by von Bahr and Esseen 
(1965, Theorem 2), there is a constant C such that 

p- 

< C2(^-")"iE 



E 



2(i-°')"i-l 

< C2P-^2(^"°^"iE 



< C2P-^2(^"")"iE 



{9niu) - g^^iu)} d 
9n{u) dn 



{9n{.u) -5(„i(n)}dn 



0,(l-ci)ni 



Ao.ri 



(1 — a)n2 



du 



(1 — Ck)?!]^ 



P1 



fi'n(M) dw 

^0,(1-q)71i 

Since /„//„i = /n-nu By the scahng property (9), we have 



0,(1 — 







Thus 

E 



gn{u) du 



0,(l~a)7i2^ 



2(l-Q)ni(V'(p<7)-p) 



E[e"'"- 




(E[e"''n 


(0)])P 






/)-p) E 


[(/" 







E 



1 \ p- 

9n~niiu) du 



gn~m (u) du 



Denote u„ = E 



Jo9n{u) du 



We have proved the following recurrence: 



Ur 



u 



n—ni ■ 



By choosing ko large enough, this yields that for any e' G (0, e), 

Un < B + 2''^~"^"''^^'^^^'^^~^'^^'^^~^'^Un-ni 

Thus, there exists a constant D such that 

This proves (37). The bound (38) follows by replacing the measure du with a discrete 
measure. □ 
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Lemma 5.3. Let < a < 1. For p > 1, there exists a positive constant C such that 

lim 2"^(^)e^('')/-^('')E[M9(Ao,„)] = m(g) , (39) 
|E[e5^'n(o)M^(Ao,„)] - E[Af«(Ao,„)]| < C2~""2-"^('') , (40) 



E 



|e'''"'"(°)M^(Ao,„) -M«(Ao,„)rl < C'2-""/22-nC(pg) 



(41) 



Proof. Define the sets /„, Bn{u), u E [0,2 "] by 



0<«<2-" 



See Figure 2 for an illustration. By definition of the function ip and the measure /i, we 
have, E[e''^(-^")] = e^^'?)'^^-^") and 



This yields 



ds dt 



^ _ 2-n /•°° 1 - 2-" , 

dt + / . . dt 



log(/„) - 2""(C' _ 1) + 1 _ 2-" = 1 - log(/„) - 2""" = /i(A,„(0)) - 2- 



E[e9^an)] = E[e''"''"(°)]{l + 0(2^"")} . 



(42) 




Figure 2: The sets J„ and Bn{u) 
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Then, wi^X'^) = P{In) + P{Bn{u)), where the two summands are independent and we can 
write 

M{Ao,n)= e"""(")M„(du) =e^('"M e^(^"("))M„(du) 



2-n 



e„ / e^(^"('^))M.(d«) , 



(43) 



with = e^*^^")M„(Ao,n) and M„(dM) = M„(dM)/M„(Ao,n) is a probabihty measure 
on Ao,„. Note that P{In), ^n(^o,n) and the integrand in the integral in (43) are mu- 
tually independent. This yields 



M'?(Ao,„,) = a I e^(^"("»M„(dn 

= a + <r Po(i?n(w))M™(dw) 
+ qa f {e^(«"(") - 1 - Po(Sn(w))}Mn(dw) 
+ e^(«"("»M.(d«)j {e^(^"("» - l}M„(d 

By elementary calculus, it is seen that for g > 1, there exists a constant Cq such that, for 
all X > —1, 



0<{l+xy -l-qx< Cg{x' + \x\''''^) 



(44) 



Applying this bound and (27) and the fact that yu(i?„(0)) = /i(i?„(2 ")) = 2 we obtain, 
for p > 1, 



E 



pq 



Sn 



gP(iJ„(-))M^(d^) - 1 - g / {e^(^'^W) - l}M„(dn) 
I Jo 



sup |e^(^"(«))-l|2p 

uG[0,2-"] 



Similarly, applying the bound (28), we obtain 



E 



{e^{i^4«))_i_p^(5^(^))}M^(dn) 



sup |e^(^"("))-l-Fo(5„(n))|^' 

«G[0,2-"] 



< C2-"" E[^P' 
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By the independence of Po{Bn{u)) and the measure M„, we also have 



E 



a I Po(i?n(w))M„(dM) 




E 



a / E[Po{Bn{u))]Mn{du) 




< 2-"-E[a] . 



We have now proved that 

E[M^(Ao,„)]=E[a{l + 0(2-"")}. 

Moreover, applying (42) and the independence properties, we have 

E[a] = E[e''^(^")]E[M;?(Ao,„)] = Eie^'^'-W] E[M^(Ao,„)]{l + 0(2'"")} . 

The last two bounds yield (40). By the scaling property, we have E[M^(Ao,n)] = 2~"^(''W(g). 
Since Efe''"""^^)] = e'^^'iHn^'^'^\ (39) follows. 

To prove (41), we use the bound (26) to obtain 



E 



cpq 



2-n 



< C2-""/2e[^p9] . 



Gathering this and the previous bounds yields (41). 



□ 



Proof of Lemma 3. 3. By the scaling property, and since the random variables -Do,fc,n,g are 
centered, we can write 



k' 



- T^'i\k - -)^(2g)cov (^M''([0, 

2 y rC 

- 2-('')A:«(2'')cov f M^([0, \\).M\\\ - i-, 1]) 



1/2])' ^^^[^-2^ 



,1]) 



+ 22-('?)(/t - l/2)'^(2g)^Q^ i ^gQo, — — -]), M^([l - 
+ 22-Wfc^(2.)cov (^^^"([0, ^]), M''([l - ^, 1])^ 
- 2^('^)A:«(2'^)cov (m^([^, ^]), M'^([1 - i, 1])) 



1 



Ik-X 



1]) 



+ 22^('^)(A;-l/2)'^(2g)^Qy /j^^gQ 



1 



"Ik - \ 1k-\ 



1 



2k - 1 



,1]) 



+ 2^^('^)fc^(2'^)cov (Af^([^, i]), Ar^([l - ^, 1])^ 
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We will prove below that each covariance terms that appear above is of order k ^'''^^^ ^, 
which yields 2'"^^'^'i^E[Do^o,n,qDo^k,n,q\ = 0(fc^(29)"2C(9)-i)^ and since C(2g) - 2C(g) = 2ip{q) - 
^{2q), the bound (13) is proved. 

Let us now prove the bounds for each covariances. The derivation is similar for all of them, 
so we will only study one explicitely, say cov([M'^([0, l/k]), M'^{[1 — 1/k, 1])). 

For k > 3, define I = 1 — 2/k and Mi{du) = e'"^^^^^ M (du) . By construction, the measure 
Ml is independent of {wi{u)} and M/([0, l/fc]) is independent of M{[1 — l/k,l]). Define 
the sets and by 

A, = Mi/k) \ Ail - i/k) , B, = Ml - i/k) \ Ai{i/k) , 

For u < k and v > 1 — l/k, define 

Ek,u = Cl'^'" \ , = A{v) \ Ail - l/k) . 
See Figure 3 for an illustration. 



k,u 



El. 



- + - 

2 ^ k 



^k 



D 



k,u 



Bu 



k,v 




Figure 3: The sets A, 5, C, D, D', E, E' 



Note that all these sets are above the horizontal line at level 1 = 1 — l/2k, hence P{A) 
is independent of Mi and P{A) is independent of P{B), where A,B are any two of these 
sets. Note that 

I I /^U,V — 1/fe 

'Ju<l/k,v>l-l/k^k ~ ^k ) 
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and let this set be simply denoted Ck- Note that Dk,u C Cfc, Ek^v C C^, Z)^ ,^ C -D^ q and 
„ C -Efe 1- We moreover have 

/i(Afc) = = 1 , 

^^{Ck) = -\og{l-2/k) , 
/xPI,o)=M^m) = 1/(A:-2). 

With this notation, we have, for u <l/k and v >1 — 1/k, 

wi{u) = P{A,) + P{Cn + ^(^L,J + P{Ek,.) , 
wi{v) = P{B,) + P{Cn + P(D,,„) + P{EiJ , 

Recall that the random measure P can be split into two independent random measures Pq 
and Pi such that P = Pq + ^i- For i = 0,1 and u G [0, 1/k], define 7ri,fc(tt) = Pi{D'^^^) + 

P,(C:''-'/') and 

7rfc(n) = 7ro,fc(n) + 7ri,fe(n) . 
Similarly, for « = 0, 1 and w G [1 - 1/A;, 1], define 7r^,(t;) = Pi(^fc_J + Pi(C;/^''') and 

Let Ml and M/ denote the normalized measures MJ Mi{\0,l/k]) and Mi/Mi{[l - 1/k, 1]). 
Denote finally 

a = M,([o,iA]), a = M([i-iA,i]) , 



Ik 



JO 

1/fc p1 



1-1/k 



Jl-l/fc 



With this notation, we obtain 

E[M^i[0,l/k])M^i[l-l/k,l])] 



E[e«^(^fc)]E[e'?^(^*)] E 



1/fc /.I 



'0 Jl-l/fc 

-1/k nl 



+ e^'^^'^^E 



'0 Jl-l/k 

Jo Jl-l/k 
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Since ttq and ttq are independent of the measure Mi, we have 



E[C^,] = E 



E[no,kiu)]Mi{du) 



< afc{/i(D^,o) + KCk)} = 0{ak/k) 



Similarly, E[^^/3fc] = 0{hk/k). Applying the bound (28) (and some algebra), we obtain 



sup sup |e-'=(")+<('') - 1 - 7ro,fc(w) - tx'^^M \ 

«G[0,l/fc] v&[l-l/k,l] 



0{akhk/k) . 



For the last term, we apply the bounds (44) and (27) and obtain 



E 



l/k rl 



I / / e-^-(^)+<('')M,(dn)MKdtO 

Jl~l/k 

.l/k pi 

1-q / {e^'=W+<(^) - l}Mi{du)Mi{dv) 

Jo Jl-l/fc 



0{akhk/k) . 



Altogether, we obtain 

E[M«([0, 1/A;])M^([1 - 1])] = e^'^^"'^ akhk{l + 0{l/k)} . 

We also have 

M([0, l/A;]) = e^(^^)a x ^ + + / {e^^(") - 1 - TXo,k{u)]Mi{du) 



(45) 



M{\l/k, 1]) = e^(^^)e. X 1 + + / {e<(^') - 1 - T,',^^{v)}M'Mv) 

V J i-i/fc 

By similar techniques, we obtain 

E[M'^([0, l/k])] = e^^'i^ak{l + 0(1/A;)} , 
E[M''([1 - l/A;, 1])] = e^^^'^il + 0(1^)} , 

and thus 

E[M^([0, l/A;])] E[M''([1 - l/k])] = e^^'^^^afc^^jl + 0(1/A;)} . 

Gathering (45) and (46) yields 

cov(Mf([0,l/J^]),Mf([l - 1/A;,1])) = 0{akhkk-^) . 

Finally, the previous bounds also imply that = 0{k~'^^''^) and bk = 0(A;~''*^''^) so that we 
have finally proved that cov(M''[0, l/A;]), M''([l - l/k, I]) = 0(A;-2C(9)-i). □ 



(46) 
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5.2 Multifractal random walk 



As a first step we require tlie following result for aj^k,n,H which is analogous to Lemma 5.3. 

Set aj,k,n,H = e"''"(*^''=)5j,fc,n,// with 

~Slk,n,H= I , I , |n-tf^-2M„(du)M„(d.;) 

k,n k,n 

and for ji ^ j2, 

Ph(Ji, J2, fc, fc ) = ^ J . 

Oji,k,n,H 0j2,k,n,H 

Lemma 5.4. For p > 1 snc/i that 2pq < and for r > 2, there exist rj,C > and 
uniformly bounded constants Cg^H{k,k') such that 

2nCH(29)g^(29);-^(2g)^[^29^^^^j _^^^2g) = 0(2""'') , (47) 

|2"C.(2.)eV^(2.);-^(2,)E[5._^_^^^5.^, ^_^] _ ,^^^^k, k')\ = 0(2-"^) , (48) 
2<"^'^Mall^^^ - all^^^l^] = 0(2-"") , (49) 

2-^"^'^MK,k,n,HaU',n,H ' <,k,n,H~<,0,k' ,n,Hn = 0(2-^) , (50) 

The proof is along the lines of the proof of Lemma 5.3 and is omitted. We can now prove 
Lemma 4.2. 

Proof of Lemma 4-2. The first part follows by Lemma 5.4, along the same lines as the 
proof of Proposition 3.1. For the second part of the Lemma, where we assume H < 3/4, 
using the notation introduced in Section 4.2, we have 

L-l 2" 

Let Qriq) r > 0, be the coefficients of the expansion of Gg over the Hermite polynomials 
{Hr,r > 0} (which are defined in such a way that E,[Hk{X)Hi{X)] = k\ if k = I and 
otherwise). Note that since G is a centered even function, gr{q) = for r = 0, L Then, by 
Mehler's formula (see e.g. Arconcs (1994)), we have 

L-22^"-(^)varM(^L,.(X,g)) = ^L^r, g) , 

r=2 ^' 

with 

(L-l 2"-l 
j=0 k=0 

^ ^-222...(.) g X;^kn(jl,J2,A:,fc')4,M,H4.^'.n,H 
ii J2=0 k,k'=0 
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for r G N, r > 2, and the conditional correlations (which are zero if if = 1/2) are 

PH.nUl, J2, k,k) = COY M[Yj^,k,n,yj2,k',n) = • 

^ji,k,n,H (^j2,k' ,n,H 

By Lemma 3.1 in Ludciia (2008), for ji < j2 and k < k', we have the bound 

PH,n{ji,J2, k, k') < min(l, C\{j2 - Ji)2" + {k' - fc)!^^"^) (51) 
for some deterministic constant C. 

We start by proving that for H < 3/4 and q < q^, there exists a constant r(r, q) such that 

lim 2"(2^(^)-^(2g)+i+x)E[r„(r, g)] = r(r, g) . (52) 

n—^oo 

By the scaling property, 

n^UnJ = 2-"^-(^)m^(g) , 

with Cniq) = qH — ip{q). Thus, denoting v^{q) = 2ip{q) — 'ipi/lq) + 1 + x, by stationarity, 
we have 

2™x(.)E[r„(r, q)] = r! m^(2g) + 2~-2<-('^^r\ ^ E[p^^^,{k, k')al ,^^^^al,,^J 

k^k' 

+ 2-(i+x)2K«(2.)H J]E[p^,,Jj,j',A;,fc')aI.,n,H4,^',n,H] (53) 
Consider the middle term. Recall that 

-(fc+l)2-" Mk'+l)2- 



PnA^. k')al,k,n,H<,k',n,H = / / l« " V^-' M (du) M (dv) 

Jk2-" Jk'2-" 



"(fc+l)2-" ^(fc+l)2-" 

X <( / / |n-t;|2^-2M(du)M(dtO 



(g-r)/2 



fc2-" Jfc2-" 
(fe'+l)2-" /•(fc'+l)2- 



X <( / / |n - v\'^"-'^M{du)M{dv) ' 

'a:'2-" Jk'2-" 

Assume that k < k' and denote £ = /c' — + 1. By the scaling property and the stationarity 
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of the increments of M, we have 

Pn,Hi^^ ^ )'^0,k,n,H'^0,k',n,H 

[Jo Ji-i/e 

.1/1 fin ^ 

'i-i/e Ji-i/e 



with 



a]= I / \u-v\''^~''M{du)M{dv) 



'0 Jo 
"1 /.I 

l2 I I l„. „.|2//-2 



Ji-i/e Ji-i/e 

aeOe 

for some deterministic constant C. Thus the middle term in (53) can be expressed as 

2n-l 

=1 

2n-l 

2^!2<«(2'?)2-"{29^^-^(2g)} ^(1 - £2-")£2'?^-'^(2^)E[Q^a^ 

271-1 

2r! XI (1 - ^2-")£'^«(2g)E[gr^g^g] 



1=1 



This sum has a hmit if the series Yl^i ^'^"^'^''^'^[Q^'A^'e] is summable. This holds true since 
applying the bound (54) and Holder's inequality yields 



1/2 



Applying the stationarity of increments and the scaling property of the MRM M yields 
E[af ] = E[6f ] = l-^H^'^i^mH{q), hence 



■2) 
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Since r > 2 and H < 3/4, the series i^^"^^ is summable, and thus 

hm y"(l - i2-^) i^"^^'^ E[Qla%l] = yf«(^'') nQ^aX] . 

t=l i=l 

Consider now the last term in (53), say i?i?„. Using the bound (51), the scahng property, 
the fact that the aj^k,n,H are 2-dependent, and H < 3/4, we have 

L 2" 

RRn < C2"^C«(2'?)-2ChW} J2 ^(j2" + kf^-^ = o(2"^2^('')-^(2g)}) _ ^(1) _ (55) 

j=i k=i 

This proves (52). 

We now prove that if -ff < 3/4, for each r > 2, 

r„(r,g)/E[r„(r,g)] ^ 1 ,a.s. (56) 

or equivalently 

2n{i+x-^(2,)+2^M}r^(^^ g) ^ r(r, g)] , a.s. (57) 
Write 2"{i+x-V'(29)+2^(9)}r„(r,g) = r!(^„,i + + ^„,i) with 

L-l 2"-l 

j=0 fc=0 
L-l 

2"-l 

0<j¥=j'<K k,k'=0 

The bound (55) and a Borel-Cantelh argument imphes that Sn.3 — )■ almost surely. The 
bound (49) implies that 

2-.(2,)^-i g - ^ a. (58) 

j=0 k=0 

Thus we consider Sn,i = 2'^'^"^^'^''L~^ ^^=0^ Sfc=o^^jlnH- Using Lemmas 5.2 and 5.4 and 
mimicking the proof of Proposition 3.1, we obtain that S'„^i — )■ mH{2q) a. s. 

By stationarity and 2-dependence in j, to deal with Sn,2) as in the proof of Proposition 3.1, 
it is enough to prove that 

PI 



E 



0<fc7^fc'<2" 



Q^2'^ex-V)n^ (59) 
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for some 77 > 0. Since all quantities involved are nonnegative, we can use the bound (51), 
and suffices to obtain a bound for 



E 



Z |Ab I %,k,n,H"'0,k' ,n, 

Q<k^k'<2" 



H 



Define 
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k- I I |M-t;p^^'M„(dM)M„(dt;) . 



Then aQ^k,n,H = Ske'^'^'"^'''^ and using the bound (50), we obtain 

^ ^ ^ l'^ '^1 l"'0,k,n,H"'0,k',n,H "o,fc,n,,H"o,fc',n,H J 

0<fc^fc'<2" 

Thus we need to obtain a bound for £[5*^4] where 



0(2' 



^">4 — ^ I ft- ft^ I "'0,k,n,H"'0,k' ,n,H ' 

0<fc7^A:'<2" 

wich we further decompose as Sn,4 = Sn,5 + 5n,6 with 

0<fc^fc'<2" 

The bound (48) implies that 2"^^^(2g)E[^9^g^]£[g<?«;,„(fc2~")+g«,,„(fc'2-")] uniformly bounded, 
thus 



2"" ^ \k-k' 

0<kj^k'<2" 



r(2H~2) ^ 



,qwi„(k2-")+qwi„ik'2-") 



]E[e9"'in('=2-")+g«)i„(fe'2-" 

Applying Lemma 5.2, for p such that 2pq < and e' < p — 1, we can prove that 



E 



2^" ^ |A;-A;' 

0<fc5>^fc'<2" 



r(2H~2) ^ 



qwi^{k2-")+qwi^{k'2-") 



]^[eqwi„ {k2-r')+qwi„ (fc'2-")j 



< (-f^-{^(2p5)-pV'(2g)-e'} ^ 



(60) 



If 2pq < g^^, we have 

(1 - a){iji2pq) - p^{2q) - e'} - ex < (1 - «)e(l + x) - e' < e - e' - ae{l + x) 
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which can be made negative by choosing e' close enough to e. 

To deal with the last term, as in the proof of Proposition 3.1 we use the conditional 2"" 
dependence of the random variables 6k- We obtain the bound 

for small some rj > 0. We have proved (59), and thus (56) holds. 
We can now define 

r=2 

As Yl'^2(^^-)~^9r{Q)'^ < 00 and r„(r, g) < r„(2,g), then by bounded convergence, the 
previous series is convergent and 

2n(2^(g)-^{2g)+i+x) ^^^^^ S L,n{X , q)) ^ T{q) , a.s. 

Hence, for 2q < q^, and 1/2 < if < 3/4, by Chebyshev's inequality and an application of 
the Borel Cantelli Lemma, we obtain (19). □ 

Proof of Lemma 4 ■ 3. For k > 1, denote 

pl/k i>l/k 

Uo,k= / \u-v\'^"~^M{du)M{dv) , 

Jo Jo 

j'l/2k j-l/2k 

U*f^= / \u-v\'^"-^M{du)M{dv) , 

Jo Jo 

nlk nl/k 

U2,k= / \u-v\^^~^M{du)M{dv) , 

Jl/2k Jl/2k 

Vo,k= [ \u-v\^''-^M{du)M{dv) , 

Jl-l/k 

V*,,= / / \u-v\^"-^M{du)M{dv) 

Jl-l/k Jl-l/k 

V2,k= [ [ \u-v\^^-^M{du)M{dv) . 

Jl-l/2k Jl-l/2k 

Then, by the scaling property, we have 

2<«(^'')E[Un,oUn,k] = k^"^''^^coviUo,k, Vo,k) - 2^"^''\k - l/2)^«(2<?)cov(f/o,,, Vi,k) 

- 2^"^'^^k^"^^'^\cOY{Uo,k, V2,k) - COv(t/i,fc, Vo^k) + COv(t/2,fc, Vo.fc)} 

+ 2^^»^''\k - l/2f"^^'i\coM{U^^k, V^i,fc) + cov(?72,fc, V^i,fc)} 

+ 2'^«^'i^k^»^^'^\cOY{U^^u. V2,k) + COv(t/2,fc, ^2,^)} . 
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All the covariance terms are of the same order, and as in the proof of Lemma 3.3, we will 
only consider the first one. Denote I = 1 — 2/k and define the measure Mi and all other 
quantities as in the proof of Lemma 3.3. By the same type of decompositions as in the 
proof of Lemma 3.3, defining 

"l/k j-l/k 



Ck,H= / / \u-v\^''-''Mi{du)Mi{dv) 
Jo Jo 



^k,H = / 1^^ - v\'''-'Mi{du)Mi{dv) , 

Jl-l/k Jl-l/k 

ak{q) = E[Cl^] , Mi{du) = a,{q)-'/^ Mi{du) , 

pl/k r-l/k 

ak= / \u-v\^^~^7rk{u) + 7rk{v)}Mi{du)Mi{dv) 



h{q) = El^j^J , M/(dn) = hiq)-'/'Mi{du) , 

^ = f f \u - ^r^-^{<(«) + 'K',{v)]M'Mu)M[{dv) , 

Jl-l/k Jl-l/fc 

it can be shown that 

mik] = + 0{k~')} + qe'^'^'^MClnS^,] , 

Wo'k] = e'^^'\{q){l + 0{k~')} + qe'^^^\{q)E[aJk] , 
mikK'k] = e'^^'^^d,{q)h{q){l + 0{k^')} + qe''^^'^^ {h{q)E[Cl^a,] + dk{q)E[aJk]} • 

This proves (21). □ 
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